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No-hair like relations between the multipole moments of the exterior gravitational field of neutron 
stars have recently been found to be approximately independent of the star’s internal structure. This 
approximate, equation-of-state universality arises after one adimensionalizes the multipole moments 
appropriately, which then begs the question of whether there are better ways to adimensionalize the 
moments to obtain stronger universality. We here investigate this question in detail by considering 
slowly-rotating neutron stars both in the non-relativistic limit and in full General Relativity. We find 
that there exist normalizations that lead to stronger equation-of-state universality in the relations 
among the moment of inertia and the quadrupole, octopole and hexadecapole moments of neutron 
stars. We determine the optimal normalization that minimizes the equation-of-state dependence in 
these relations. The results found here may have applications in the modeling of X-ray pulses and 
atomic line profiles from millisecond pulsars with NICER and LOFT. 
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I. INTRODUCTION 

One of the most interesting physical results that one 
may derive from neutron star (NS) observations is a bet¬ 
ter understanding of the supra-nuclear equation of state 
(EoS), i.e. the relation between internal density and pres¬ 
sure at densities beyond nuclear m- Imagine, for ex¬ 
ample, that one were to observe the NS mass and its ra¬ 
dius independently. Since the relation between the mass 
and radius of a NS is highly sensitive to the EoS, such 
observations would place strong constraints on the lat¬ 
ter. But the NS radius, in particular, is currently very 
hard to measure with sufficient accuracy, which makes 
constraints on the EoSs not sufficiently strong [HE]. 

Ignorance of the EoS can have a strong impact on the 
amount of information that can be extracted from astro- 
physical observations. When the EoS is unknown, more 
model parameters are typically needed to fit and inter¬ 
pret the data. For example, in gravitational wave (GW) 
astrophysics, waveform models are constructed to match- 
filter the data. Such models are characterized by a set of 
physical parameters that describe the system that gen¬ 
erated the GWs, e.g. for NS binaries, these include the 
masses, the spin angular momenta, the quadrupole mo¬ 
ments and the tidal deformabilities. If the EoS is known, 
one does not need to fit for the quadrupole moments 
and the tidal deformabilities, since these are functions of 
the masses and spins only. Lacking precise knowledge of 
the EoS, these quantities must be included in the model 
parameter list and fitted for, which then dilutes the ac¬ 
curacy to which other parameters can be measured. 

The aforementioned problem can be alleviated if EoS- 
independent relations between the moment of inertia (/), 
the tidal Love number and the quadrupole moment (M 2 
or Q) of NSs can be found: enter the I-Love-Q rela¬ 
tions [i [7]. Although these relations are not exactly 
EoS-independent, they are approximately so, with vari¬ 
ations only at the percent level. After their initial dis¬ 


covery, the I-Love-Q relations were confirmed using dif¬ 
ferent EoSs [5], using binary systems [S], using magne¬ 
tized NSs [To] and proto-NSs m , allowing for rapid rota¬ 
tion [HHHj, through a post-Minkowskian expansion HZ], 
using NSs with anisotropic pressure |I8j and in alterna¬ 
tive theories of gravity [ilTlIIMll- An important ex¬ 
tension of the I-Love-Q relations was also recently found 
that relates all of the multipole moments of the exterior 
gravitational field of NSs to just the first three in an ap¬ 
proximately EoS independent fashion Oils]. 

The I-Q relations may help in the independent deter¬ 
mination of the NS mass and radius [24U28] from obser¬ 
vations of the X-ray pulse profiles of millisecond pulsars 
with NICER [Zni and LOFT [301 El] • In principle, the X- 
ray profile depends on the independent parameters (M, 
R, /, M 2 , f, (.), where M is the mass, R is the radius, / is 
the spin frequency and t is the inclination angle. The I-Q 
relation can be used to eliminate M 2 in favor of /, and for 
stars with compactnesses C = M/R G (0.1, 0.4), the I-C 
relation is also approximately EoS independent. Such a 
reduction in the model parameter space may allow one 
to extract the mass and radius independently from X-ray 
observations with perhaps a 5% accuracy |2ll ES] . 

Given these accuracy requirements, one must ensure 
that systematic errors are under control. One such source 
of error is in the EoS variability of the TQ relations at 
the percent level. Another is in the effect of higher-order 
multipole moments, like the octopole moment S 3 and the 
hexadecapole moment M 4 . These moments do not need 
to be included as new parameters in the pulse profile 
model because of the S3-M2 and the M 4 -M 2 relations 
discovered in [111123], which are also approximately EoS 
universal. The EoS universality in the latter, however, 
is weaker than in the I-Q relations, with EoS variability 
at the 15% level, which may introduce systematic errors 
that could contaminate the extraction of the mass and 
radius m- 

The EoS variability of the TQ relations and the rela- 
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tions between higher multipole moments depends sensi¬ 
tively on how they are adimensionalized. This is related 
to the point made by [T^] and [T3], who discovered that 
the TQ relations are approximately universal for rapidly 
rotating NSs provided one fixes a dimensionless spin pa¬ 
rameter as one explores a sequence of NSs. If one fixes the 
dimensional spin frequency, the EoS variability greatly 
increases, as discovered in [32]. This, of course, does not 
mean that the EoS universality is lost, but rather that an 
inappropriate spin parameter was fixed in the sequence 
of stars. Such findings provide clear evidence that the 
EoS variability is sensitive to the parameters chosen to 
construct the relations. 

Can one then find a better set of normalization con¬ 
stants that further reduces the EoS variability in the TQ, 
M4-M2 and S3-M2 relations so that they can be used in 
the modeling of the X-ray pulse profile? The answer to 
this question is yes and it is the main topic of this pa¬ 
per. When the I-Love-Q relations were discovered |7|, the 
moment of inertia was adimensionalized via I = 
the quadrupole moment via M 2 = —M 2 /{M^~x^), the oc- 
topole moment via = —Sz/{M'^^) and the hexade- 
capole moment via M 4 = M 4 /where x = 
is a dimensionless spin-parameter. Si = in is the spin 
angular momentum, and n is the spin angular frequency. 
With such choices, the I-Q, M4-M2 and S3-M2 rela¬ 
tions have an EoS variability of 2%, 15% and 8% re¬ 
spectively m for slowly- or rapidly-rotating NSs and 
realistic EoSs [551l41| . If instead of this choice, one nor¬ 
malizes the multipole moments via 
M 2 "®"' = -M 2 /{M^x^C), = Ss/iM^x^C^'l and 

then the maximum EoS vari¬ 
ability in the new I-Q, M 4 -M 2 and S 3 -M 2 decreases by 
a factor of 2 or 3, down to ~ 1%, ~ 6%, and ^ 2% 
respectively, for slowly-rotating NSs with the same EoSs. 

The above choices of normalization are not unique, and 
in fact, there is an entire family of normalizations that 
minimizes the degree of EoS variability. To find this fam¬ 
ily, we divide our study in two parts: (i) an analytic, non- 
relativistic treatment, where we consider polytropic EoS 
with index n S (0,1); and (ii) a numerical, fully relativis¬ 
tic analysis where we consider realistic EoS [53H3T| . In 
both cases, we focus on slowly-rotating stars with masses 
M G (I,2.5)M0 in the Hartle-Thorne formalism [42l 143] 
and without magnetic fields, as this is appropriate for 
recycled pulsars, even with millisecond spin periods. We 
compute the multipole moments for these stars and then 
adimensionalize them in the same way as originally dis¬ 
covered in [7] and explained above, but with an extra 
factor of C'“*’, where the power Op is different for each mo¬ 
ment. Each of the TQ, M 4 -M 2 and S 3 -M 2 relations then 
depends on two free parameters, (a/,aM, 2 ), 
and (as_ 3 ,aM, 2 )- We discretize this space and compute 
the maximum EoS variability at each point to find the 
set that minimizes it. In both cases, this set can be de¬ 
scribed by a straight line in {aj, qm, 2 ), 4 ,^ 114 , 2 ) and 

(as, 3 iaM, 2 ) space, as shown in Fig. fd The slope and 
y-intercept of this line is slightly diiFerent in the non- 


relativistic and in the relativistic calculations. Observe 
that the best normalization does not always agree with 
that originally chosen in [7] . Similar results are found for 
the S 3 -M 2 relation. 

The remainder of this paper deals with the details of 
the results discussed above and it is organized as follows. 
In Sec. [n] we review the original TQ and three-hair re¬ 


lations for NSs, as derived in [T3]. In Sec. Ill 


we re¬ 


formulate the equations in the non-relativistic limit with 
a generic normalization to study how the EoS universal¬ 
ity depends on this. In Sec. EYl we repeat this analysis 
but in full General Relativity. In Sec. [^ we conclude 
and highlight possible directions for future research. All 
throughout, we use geometric units in which G = 1 = c. 


II. ORIGINAL THREE-HAIR RELATIONS 


In this section, we briefly review the three-hair rela¬ 
tions for NSs derived in US). The exterior gravitational 
field of an isolated and stationary mass distribution can 
be written in terms of a multipole moment decomposi¬ 
tion |44] . In the non-relativistic limit (i.e. in the per¬ 
turbative weak field, when we expand all expressions to 
leading order in powers of compactness) and in a slow 
rotation expansion (i.e. an expansion in powers of the 
product of the mass and the spin angular frequency), the 
leading-order expressions for the mass and mass-current 
multipole moments are |45j 


pir pRtie) 

Mf = 2 tt / / p{r, 6 )Pi{cos 6 ) si \4 9dO dr 

Jo Jo 


( 1 ) 


and 


Si = 




i + i 


p{r, 6 ) 


'0 Jo 


dcos9 


sin 9d6 r ^ dr , 

( 2 ) 


where p{r,9) is the density, i?*(0) is the stellar surface 
profile, n is the spin angular velocity and Pi{cos9) are 
the Legendre polynomials of order £. These moments are 
distinct but algebraically related to the Geroch-Hansen 
moments [461 - I48] and the Thorne moments [49] in the 
non-relativistic limit misD]. 

In the non-relativistic limit, slowly-rotating NSs can 
be well modeled in the elliptical isodensity approxima¬ 
tion of |ST]. In this scheme, one assumes that isodensity 
surfaces are self-similar ellipsoids and that the density 
profile is the same as that of a non-rotating star of the 
same volume as the rotating star. Within this approxi¬ 
mation, and in a suitable coordinate system, the angular 
and radial integrals of Eqs. Q and ([^ can be separated 
in terms of certain angular integrals leji and certain ra¬ 
dial integrals Rg, namely 

Mg = 2TrIg^3Rg , (3) 
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FIG. 1 . (Color Online) EoS variability in the I-Q (left) and M4-M2 relations (right) for different values of the normalization 
constants {aj,aM,2) (left) and {aM,4,CLM,2) (right). The moment of inertia is normalized via / = I the quadrupole 
moment via M2 = —M2/(M®x^C'““'2), and the hexadecapole moment via M4 = M4/(M®x"'C““’'‘)- For a given choice of 
normalization, the contours show the maximum percent EoS variability in the I-Q and M4-M2 relations computed with variety 
of realistic EoSs | 33 H 41 |. relative to a representative realistic EoS [^. For comparison, we also present the set of (a/, aM,2) and 
{aM,i,CLM,2) that leads to the least EoS variability in a non-relativistic analysis (solid yellow line) for polytropic EoSs, as well 
as the original normalizations chosen in [B] [T] (TT] (white solid circle). Observe that even in the fully relativistic calculation, 
the set of (a/,OM.2) and {aMA,o-M,2) that minimizes the EoS variability is a straight line in this space, where the slope and 
j/-intercepts are different from that obtained with a Newtonian analysis. Observe also that there are choices of normalization 
that lead to less EoS variability than the original normalization of lilZlITl]. 


and 

AttE 

Si = _(_ ^^(Ft- 1,5 - li+i,3)Ri+i ■ (4) 

The angular integrals can be calculated in closed form 
and are [5^ 

= , (5) 

and 

T T -( A 

^(^ + 2 ) ^ 1 -" " • 

, , ( 6 ) 

The radial integrals can be rewritten using the Lane- 
Emden function [53] as 

Re = Pc(^) Rn/ , (7) 

where 9 = (p/pc)^^" is the Lane-Embden function, pc 
is the central density and n is a constant, later to be 
identified with the polytropic index. The reduced radial 
integral is defined as 

nn,i= , ( 8 ) 

Jo 


where M{= Mq) is the stellar mass of the non-rotating 
configuration, ^ = {^i/ai)f is the dimensionless radial 
coordinate, with oi the semi-major axis of the ellipse, and 
^ is the dimensionless stellar radius. Observe that 
all of the EoS dependence is here encoded in the r2(e) 
function and the radial integrals Rg. Notice also that 
all the results reviewed above assume a single polytropic 
EoS, i.e. an EoS of the form p = but they can 

be easily extended to piecewise polytropes [TB], which 
accurately approximate realistic EoSs [H [M] . 

With all this information at hand, the mass and mass- 
current moments can be written as 


M2i+2 = 


(^_lY+1^2i+2M2i+0Tl^2i+2 
(2£ + 3) (l-e2)^ c^t+2 


(9) 


and 

^ (-1)^ 2 Q(e) e2^ M^^+^TZn, 2 e +2 

*^2^+1 ^+1 no , A ’ 

(2€ + 3) (l-e2)— ^2i+4 ^2e+2 |^/(^^)| 

( 10 ) 

where C = M/R is the stellar compactness. In the ellip¬ 
tical isodensity approximation, the angular frequency is 
simply m 


m = 2 ^? 




-I 1/2 


(5 - n)M2 7^„,2 


/(e) 


( 11 ) 
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where /(e) = [—6e ^(1 — e^) + 2e ^(1 — — 

2e^) arcsin(e)]^/^. 


Let us now derive the three-hair relations. First, we 
adimensionalize Eqs. Q and ( [To| ) via 




Ml 






Si 




, (12) 


where x = Si/M'^ is the dimensionless spin parameter. 
Then, we set £ = 0 in the adimensionalized M 2 i +2 equa¬ 
tion to solve for compactness C as a function of M 2 - And 
finally, we substitute this expression back both into the 
adimensionalized version of Eq. and Eq. ( |Io| ) to obtain 
M 2 i +2 and § 21+1 as functions of M 2 - These expressions 
can then be expressed via the single equation 


M2i+2 + iS2i+i — Bn,eM2{M2 + iSi) , (13) 


where 


Bn,i — 


(2£ + 3)1A <+2'^X.'//+2 


31+1/^ 


l^'(ei)l?? 


(14) 


This is the so-called three-hair relation for NSs in the 
non-relativistic limit, which is approximately EoS inde¬ 
pendent. Indeed, all of the EoS dependence is encoded 
in the coefficients, which were shown to vary weakly 
with n for any given £ |15j . Similar expressions hold in 
the relativistic regime for rapidly rotating stars [HE]. 

The I-Q relations can be derived from the expressions 
presented above. Setting £ = 0 in Eqs. (|^-([TT|), we then 
obtain 


M 2 = 


le^ 

V 


(15) 


This equation, however, is not yet the I-Q relation we 
want, because the right-hand side still depends on e and 
y. We can solve for y as a function of e by setting £ = 0 
in Eq. (10), solving for compactness and using the result 


in the definition of y = S'l/M^. Doing so, we find 


X = 


715 


Jl/4 


4 (1- 671 / 44 /^ 


/(e), 


(16) 


and substituting this into Eq. (151 we obtain the I-Q 


relation. Expanding this in the slow-rotation limit, we 
can write the I-Q relation as 


M 2 = \riA 


n,0 5 


(17) 


where 


An ,0 — 


25(5-n) 
1152 


21 


27 TZl 




l^'(7)|- 


(18) 


III. TOWARD EXACT UNIVERSALITY: 
ANALYTICAL RESULTS IN THE 
NON-RELATIVISTIC LIMIT 


Can the approximate EoS universality of the NS no¬ 
hair relations be improved by modifying the normaliza¬ 
tion of the multipole moments? Let us then re-define the 
dimensionless multipole moments as 


and 


M, 2 ^ + 2 .) _ {-lY+^M 2 i +2 

2^+2 J^2l+‘i^2i+2(JaM,2e+2 


o(as,2t+l) _ _ ( _4)_ ^ _ 

‘^21+1 — j^2l+2.^2l+l(jas,2t+i 


(19) 


( 20 ) 


where aM, 2 f-i -2 and as, 21+1 are real constants that provide 
flexibility in the adimensionalization. Notice that to re¬ 
cover the original normalization, one must set aM, 2 i +2 = 
0 and as, 2 i.+i = 0- 

The approximate no-hair relations can now be ob¬ 
tained with the new normalization. First, we eliminate 
the compactness from the dimensionless multipole mo¬ 
ments by solving Eq. (19) with £ = 0 for C in terms 
of Substituting this back into Eq. (19) and 

Eq. (20), we find 


Xr{aM,2l+2) _ 4°^ „_a+f+l f Mn \ ^ „2\-a/d, 

^ki +2 -2£ + 3 [ej ^ ^ 

X |ll'(6)|'-“ Kr2'-' ^n.2£+2 

X , (21) 

c(“S,2/!+l) _ 4^ f MM\ ^ _ 2',-/3/3 

*2^+1 “ 2£ + 3 V e j ^ ^ 

X |7(a)r-^ K:2~" T^n,2i+2 

X , (22) 


where a and /3 are new constants defined by 

^ _ (£ -I- l)aM,2 — aM,2£-t-2 
—aM ,2 + 2 
p _ £aM,2 H—as,2^-i-i 
~«M ,2 + 2 


(23) 

(24) 


Unlike in the original no-hair relations, summarized in 
Sec. [TT| this time the multipole moments depend on the 
spin angular frequency Q and the eccentricity e. One can 
eliminate this dependence by choosing om, 21-12 = (£-!- 
l)aM,2 and as,21+1 = £aM,2, which leads to a = 0 = / 3 , 
and thus to 


rAaM. 2 e+ 2 ) _ 3^+1 |i?'(gi)|^ ef T^n,2i+2 


M 


2 e-i 2 


(2£ + 3) <+' 


and 


Q(as,2i+l) 

*^2^+1 


3^+1 |t^'(/i)l^ eY n^,2i+2 

(2£73)<+i 


(25) 


(26) 
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These two expressions are similar to Eqs. (10) and (11) 
of [15], but here, they are obtained for a larger class of 
normalizations, that in particular, includes the original 
normalization of m- 

With these general expressions at hand, let us now in¬ 
vestigate a few examples of universality, beginning with 
the I-Q relations. Generalizing Eq. (22) of (TS], we adi- 
mensionalize the moment of inertia via 


jiai) _ 


Si 




(27) 


where a/ is a new normalization constant. Proceeding in 
the exact same way as in the previous section, the ex¬ 
pression for X, still defined via y = as a function 

of eccentricity is 


X = 


•\/3(l — e^) 9 (5 _ 

2-i+6a:5i-2x^-|+2^ 


/(e) 


(28) 

where x = g^^^°)^ and An/ is still given by Eq. (181. 
Clearly, we recover Eq. (16) when a/ = 0, and thus x = 
5/12. Solving for compactness from the £ = 0 version of 
Eq. ( plj ), and using Eq. (28), we get the 1-Q relation in 
the slow rotation limit: 


MAm.2 ) 


2 \ ^“ 1 +^ 
3, 


ai-2aM,2 

'TP 2aj+4 

'^n,2 


i^'(a)r 




7+4 


Cl 


An,0 i 


Notice that the new I-Q relation depends only on the 
set (qm/tCLi). Clearly, we recover Eq. 0 when gm/ = 
ai/2, which includes the original normalization of [15] . 

Given the new I-Q relation in Eq. ( |29| ), we can now 
address whether there are particular sets of parameters 
{aM. 2 ,ai) for which the EoS universality is improved, 
i.e. for which the variability of the I-Q relation with dif¬ 
ferent EoSs is decreased. One way to assess this variabil¬ 
ity is to compute the relative fractional difference in the 
I-Q relation between an n = 0 and an n = I polytropic 
EoS. This is analytically tractable because exact solu¬ 
tions to the equations of stellar structure exist for such 
polytropic EoSs. This fractional difference reduces to 


M, 


(fiM,; 


|n=0 ~ M, 




\n=l 


M. 


(aM,2) I 


1 - 4 • 5” 


2 |n=0 

“M,2+“/+3 ^{^M,2 + ^') 


1+2 


TT 


1+2 


TT^ — 6 


“M,2 + l 

a ,+2 


(30) 


and setting it to zero, we can determine the parameters 
{aM, 2 ,CLi) that lead to exact universality between these 
two EoSs. 

Eigure presents contours of the relative fractional 
difference, including the zero contour, in the {aM, 2 ,ai) 
plane. Observe that there exists a one-parameter family 
of {gm/tCIi) for which the relative fractional difference 
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FIG. 2. (Color Online) Relative fractional difference in the 
Newtonian I-Q relations for different {aM, 2 ,cti) normaliza¬ 
tions. The fractional difference is here calculated by comput¬ 
ing the TQ relations with various polytropic EoSs relative to 
a polytropic EoS with index n = 0.643. The set of {gm, 2 , 0 - 1 ) 
in the violet region leads to the most EoS independence in 
the I-Q relations, while the set in the red region leads to a 
variation of the relations with EoS of ~ 50%. The green line 
shows the choices of {aM, 2 , oj) for which the relative fractional 
difference is not well-defined. The white regions lead to EoS 
variability in the TQ relations that exceeds 50%. The yellow 
line shows the one-parameter family of (aAr, 2 ,ar) that leads 
to zero relative fractional difference in the TQ relations when 
computed with an n = 0 and an n = 1 polytropic EoS. The 
original normalization of laizi is shown with a white circle. 
Observe that this normalization lies very close to the optimal 
one-parameter family. 


is exactly zero. Observe also that the original normal¬ 
ization is very close to this one-parameter family. A 
third interesting point is that the universality deterio¬ 
rates rapidly for values of the normalization constants 
away from the one-parameter family. Finally, observe 
that there are values of (aM, 2 ,a/) for which Eq. (30) di¬ 
verges, e.g. if you take the limit as a/ —f —2 from the left, 
Eq. (30) diverges to positive infinity. This is because in 
such a limit, does not depend on C, which in turn 
means that one cannot solve C for which is crucial 

in deriving the TQ relations. 

These results suggest that one may find some choices 
of normalizations for which the TQ relations are perfectly 
universal with respect to any EoS. To investigate this fur¬ 
ther, let us re-compute the relative fractional difference 


^ Although not easily seen in the figure, the original normalization 
does not actually lie on the line. 






























where (5i = and 

^ aM,2 + l’ 
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FIG. 3. (Color Online) Values of (aAf, 2 ,n/) that lead to zero 
relative fractional difference in the I-Q relations calculated 
with polytropic EoSs with n € [0,1] and a reference poly¬ 
tropic EoS with n = 0.643. The quadrupole moment and 
the moment of inertia are computed by numerically solv¬ 
ing the Lane-Emden equation. Observe that all lines cluster 
around each other. Observe also that the original normaliza¬ 
tion [(aM, 2 ,a/) = (0,0)] used in [SIE] lies very close to any of 
the lines. 


with polytropic EoSs of different indices. Given any poly¬ 
tropic index, we numerically solve the Lane-Emden equa¬ 
tion and compute the relative fractional difference of the 
I-Q relations with that index and a reference EoS. For the 
latter, we choose a polytropic EoS with index n = 0.643, 
which is the average index of all of those presented in 
Table I of [TB]; the latter accurately approximate realis¬ 
tic EoSs that support NSs with masses greater than 
2Mq [56]. Figure shows the one-parameter family of 
{o-M, 2 iO‘i) for which this relative fractional difference is 
exactly zero. Observe that all the one-parameter fami¬ 
lies are clustered around a single line. This indicates that 
there exists an optimal set of (cm, 2 , o/) for which the I-Q 
relation is almost exactly independent of the EoS. Ob¬ 
serve also that this optimal set happens to be very close 
to the (aM, 2 , o/) = (0, 0) point, which corresponds to the 
original normalization of the I-Q relations [Tllll5l|23|. 

Let us now investigate another example of EoS uni¬ 
versality as a function of the normalization constants, 
focusing on higher multipole moments. In particular, let 
us focus on the M 4 — M 2 and <§3 — M 2 relations, which 
can be rewritten as 


|d'(6)IC?7^n.4M2^ 


(31) 


' ^ ' ’ 

'^n,2 


(32) 


where S 2 = using Eqs. (pTj) and (|22|) with i = 1 . 


As before, we compute the relative fractional difference 
in these relations between an n = 0 and an n = 1 poly¬ 
tropic EoS, since exact solutions to the stellar structure 
equations exist in these cases: 




,-25i 


625 (tt^-6)' 


X TT^ (625 X 22'*! - 336 X 5'^4) 

- (625 X - 448 x 5^^+^) 

-b 9 ^625 X - 896 x 5*^4+^) 


and 


cifas.a) I 


s. 


|n=0 


- s. 


q(°-S,3} I 


2-2^2 


c(“S,3) I 


s: 


125 (7r2 - 6)" 


X TT^ (125 X 2^'^" - 84 X 5^") 

-b TT^ (336 X 5^=^+^ - 375 X 

-b 9 (125 X - 224 X 5'^"+^) 


(33) 


(34) 


Just as before, we now find the values of (aM, 4 , o,m, 2 ) and 
(ns, 3 ,nM, 2 ) which minimize the degree of EoS variabil¬ 
ity. Figure shows contours of fixed relative fractional 
difference in the (aM, 4 , a,M, 2 ) (left panel) and {as, 3 , aM. 2 ) 
(right panel) planes. Observe again that there are one- 
parameter families for which the fractional difference is 
exactly zero. This time, however, the original normaliza¬ 
tions of m are not quite on these one parameter fami¬ 
lies. This is important because it indicates that there are 
better choices of normalization that minimize the EoS 
variability further. Also as before, observe that the EoS 
variability increases as one chooses normalization values 
away from these one-parameter families. Moreover, there 
is a line in the (aM, 4 ,aM, 2 ) and ( 05 ^ 3 , um, 2 ) planes for 
which the fractional differences are not well-defined. 

As in the I-Q case, one may wonder whether these one- 
parameter families change if one computes the relative 
fractional differences in the M 4 -M 2 and S 3 -M 2 between 
other EoSs. Let us then repeat the study carried out in 
the I-Q case with polytropic EoSs. Figure shows the 
one-parameter families for which the relative fractional 
difference in the M 4 -M 2 and S 3 -M 2 relations is exactly 
zero between a polytropic EoS with index n and a ref¬ 
erence polytropic EoS with n = 0.643. Observe that all 
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FIG. 4 . (Color Online) Same as Fig.j^but using the M4-M2 in the (om,2, aM, 4 ,) plane (left) and S 3 -M 2 in the (as, 3 , aM,2) plane 
(right) relations. Observe that the original normalizations of [6l [T] [ 14 ] are close but do not quite coincide with the choices of 
normalization constant that minimizes the EoS variability in the relations. 


the one-parameter families cluster around single curves. 
Moreover, the original normalization chosen in m does 
not lie on any of these families. This suggests that there 
may be an optimal one-parameter family of normaliza¬ 
tions that is different from the original choice of [T3] 
which improves the EoS universality of the M 4 -M 2 and 
S 3 -M 2 relations. 


IV. EOS UNIVERSALITY 
IN FULL GENERAL RELATIVITY 


We have seen that it is possible to construct TQ, M 4 - 
M 2 and S 3 -M 2 relations that are properly normalized 
such that they are essentially EoS insensitive in the New¬ 
tonian limit and for simple polytropic EoSs; but does 
such a choice of normalization also lead to EoS insen¬ 
sitivity in the relativistic regime and for realistic EoSs? 
This section is devoted to answering this question. 

Let us start by reviewing how the TQ, M 4 -M 2 and 
S 3 -M 2 relations are computed in full GR, following [T] 
[HET]. We construct solutions to the Einstein equations 
that represent isolated, unmagnetized and slowly rotat¬ 
ing NSs, perturbatively in the ratio of the spin angular 
momentum to the stellar mass squared with the Hartle- 
Thorne framework |42[ I4d| . In order to compute the hex- 
adecapole moment M 4 , we must retain terms up to fourth 
order in the small-rotation parameter. The matter sec¬ 
tor is modeled as a perfect fluid with realistic EoSs, such 
as APR [33|, SLY |M1|35|, LS220 [36l[32|, Shen [38l|39|, 
WFFl [40] and ALF2 |4T]. All these EoSs support NSs 
with masses greater than 2Mq, which is needed given 
the two recently discovered massive pulsars (SHI EH] ■ We 


also include a few simulations with polytropic EoSs for 
comparison with the previous section. 

The equations of structure are then solved order by or¬ 
der in the slow-rotation expansion (see e.g. [71II1E7]). At 
any given order in rotation, one must solve the structure 
equations numerically in the interior of the star and then 
match them at the stellar surface to an exterior solution. 
The numerical calculations are done with an adaptive 4‘^ 
order Runge-Kutta integrator [SH]. The boundary con¬ 
ditions at the stellar center are obtained through a lo¬ 
cal analysis of the structure equations at the stellar core, 
while the boundary conditions at spatial infinity are fixed 
via asymptotic flatness. The matching of the solutions at 
IVth order in rotation yields the Vth multipole moment 
of the NSs. 


We now investigate the no-hair relations in the rela¬ 
tivistic regime and with realistic EoSs, using the numer¬ 
ical framework described above. Let us first focus on the 
I-Q relations and let us normalize the moment of inertia 
and the quadrupole moment as in Eqs. (27) and (19). 
These relations then become a function of (aM, 2 ) 0 , 1 ) and 
we wish to determine the set that minimizes the degree 
of EoS variability. To do so, we will compute the relative 
fractional difference in the TQ relations with different 
realistic EoSs, with respect to the relations computed 
with a reference EoS, which in this case we take to be 
LS220. Figurej^shows the contours of maximum relative 
fractional difference in the TQ relations in the {aM, 2 , 0 - 1 ) 
plane, i.e. for a discrete set of values in the {aM, 2 ,oi) 
plane, we compute the relative fractional differences in 
the TQ relations calculated with an LS220 EoS and all 
other realistic EoSs, and then, we create contours of the 
maximum values of those relative fractional differences. 
Observe that there still exists a set of normalizations for 
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FIG. 5. (Color Online) Same as Fig.j^but using the M 4 -M 2 in the (aM, 2 , am, 4 ) plane (left) and S 3 -M 2 in the (as.s, cim, 2 ) plane 
(right) relations. Observe again that all the curves cluster around each other and that the original normalizations of |14| do 
not quite coincide with the normalizations that lead to the least EoS variability. 



-4 -2 0 2 4 

ai 


FIG. 6. (Color Online) Same as left panel of Fig.[l|but zoomed 
out. Observe that the original normalization of |3[Z] is in the 
region of {aM, 2 ,ai) that lead to the least EoS-variability. Ob¬ 
serve also that the one-parameter family of normalizations 
that lead to the least EoS variability in a Newtonian treat¬ 
ment (yellow line) does not coincide with the regions that 
lead to the least EoS variability in full GR (purple regions). 
Finally, observe also that the region that leads to the most 
EoS variability (white region) is shifted from what one would 
expect from a Newtonian analysis (green line). 


which the maximum EoS variation is ~ 1%, as shown 
by the violet regions of the figure. Observe also that the 
original normalization used in [7], denoted by a white 
dot in Fig. is close to the best choice of normalization. 
The one-parameter family of normalization that gives the 
least EoS variability in the Newtonian case (yellow line) 


disagrees with that which minimizes the EoS variability 
in the relativistic case. Yet still, in the relativistic case, 
it seems like there is a one-parameter family that mini¬ 
mizes the EoS variability, except that now the relativistic 
corrections modify the slope of the Newtonian relation. 
Another important observation is that there exist choices 
of normalization (white region) for which the maximum 
EoS variation is not well behaved numerically; this agrees 
with the divergent region (green line) in the Newtonian 
case. A zoomed-in version of the first quadrant of this 
figure is shown in Fig. (left panel). In conclusion, we 
found a family of normalizations that lead to strong EoS 
insensitivity in the I-Q relations, which includes the orig¬ 
inal normalization of Ena, in the relativistic regime with 
realistic EoSs. 

Let us now investigate higher multipole order, no-hair 
relations and attempt to determine the normalizations 
that lead to the strongest EoS universality. In particu¬ 
lar, let us focus on the M 4 -M 2 and S 3 -M 2 relations in 
the relativistic regime and with realistic EoSs. Figures 
(right panel) and show the contours of maximum rel¬ 
ative fractional difference in these relations due to EoS 
variation in the plane of the normalization parameters. 
Observe that in both cases there is a region in the normal¬ 
ization plane that lead to the least EoS variability and 
that resembles a one-parameter family. Observe also that 
this one-parameter family in the relativistic case is dif¬ 
ferent from the one found in the Newtonian case (yellow 
lines), just as in the I-Q relations. This time, however, 
the Newtonian one-parameter family seems to be close 
to its relativistic version. Importantly, observe that the 
original normalizations in [Til [T5] are not that close to 
the relativistic one-parameter family that minimizes the 
EoS variability. For example, the best choice of normal¬ 
ization for the M 4 —M 2 relations leads to a maximum EoS 
variation of ~ 6% El: while the maximum EoS variation 
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<^ 5,3 

FIG. 7. (Color Online) Same as right panel of Fig. but for 
the S3-M2 relation in the ( 03 , 3 , am, 2 ) plane. Observe again 
that the one-parameter family of normalizations that mini¬ 
mizes the EoS variability in the Newtonian regime (yellow 
line) does not coincide with the set of normalizations that 
does so in full GR (purple regions). Importantly, observe that 
the original normalization chosen in dl is quite far from the 
purple region. This implies that there are better choices of 
normalization that would make the S3-M2 relations more EoS 
universal. 


using the original normalization was ~ 9%. Similarly, the 
best choices of normalization for the S 3 — M 2 relations 
lead to a maximum EoS variation of ^ 2%, while the orig¬ 
inal parameterization lead to a variation of 3.5% [14) . In 
conclusion, we have found the set of normalization pa¬ 
rameters that minimizes the EoS variability, and this set 
is different from the original normalization used in |I4j . 
enhancing universality by a factor of 2-4. 


V. DISCUSSION 

The recently-discovered, approximately EoS insensi¬ 
tive, TLove-Q relations are important for future obser¬ 
vations. In particular, the TQ relations are helpful in 
determining the mass and radius of NSs from the X-ray 
pulse profile emitted by millisecond pulsars. The NICER 
and LOFT collaborations hope to measure the mass and 
radius to 5% accuracy, and thus, one requires TQ rela¬ 
tions that are EoS-insensitive to at least this level. More¬ 
over, the expected pulse profile may also be affected by 


higher-order-in-spin corrections, for example induced by 
the octopole and hexadecapole moments of the NS, pro¬ 
vided the latter is spinning sufficiently fast [2]. Thus, 
one would wish to also find M 4 -M 2 and S 3 -M 2 relations 
that are the least EoS sensitive as possible. 

In this paper, we have studied how to construct TQ 
and higher multipole order, no-hair relations for NSs that 
maximize the degree of EoS insensitivity. To do so, we 
investigated whether one can normalize these quantities 
in a way that strengthens the EoS insensitivity relative 
to the original normalizations used in [7]. We have found 
the optimal one-parameter family of normalizations that 
leads to the strongest EoS-insensitivity both with a New¬ 
tonian and a fully relativistic analysis. We have observed 
that this one parameter family is strongly affected by rel¬ 
ativistic corrections, and the best normalization in the 
Newtonian limit does not necessarily lead to the least 
EoS variation in the relativistic regime. Yet,regarding 
the TQ relations, there is little gain in EoS insensitivity 
when using these normalizations relative to the original 
ones of [7] , because the latter happened to coincidentally 
lie very close to the optimal family. On the other hand, 
relativistic corrections to the M4-M2 and S3-M2 relations 
greatly affect this one-parameter family; in fact, there are 
choices that can improve the degree of EoS-insensitivity 
relative to the original normalizations of [ 2 ]. 

One could extend the present study to other modi¬ 
fied theories of gravity, such as dynamical Chern-Simons 
gravity (DCS) [6OI El]. Modified theories typically in¬ 
troduce dimensional coupling constants that characterize 
the degree to which they differ from GR. One could thus 
study what normalization minimizes the degree of EoS 
variability in such modified theories. 

Quite recently the I-Q relations were studied for highly 
magnetized NSs nnj. The authors showed that for strong 
magnetic fields {B > lO^^G) with a twisted-torus config¬ 
uration the relations lose universality. The NSs consid¬ 
ered in m are characterized by dimensional magnetic 
fields but we know that the universality of the TQ re¬ 
lation depends heavily on how we normalize the NS ob¬ 
servables. Future work could concentrate on investigat¬ 
ing whether the universality holds for highly magnetized 
NSs if a different choice of normalization for the magnetic 
field strengths is made. 
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